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Abstract. A long Gamma-Ray Burst (GRB) was detected with the instruments on board the INTEGRAL satellite 
on January 31 2003. Although most of the GRB, which lasted ~150 seconds, occurred during a satellite slew, the 
automatic software of the INTEGRAL Burst Alert System was able to detect it in near-real time. Here we report 
the results obtained with the IBIS instrument, which detected GRB 030131 in the 15 keV - 200 keV energy range, 
and ESO/VLT observations of its optical transient. The burst displays a complex time profile with numerous 
peaks. The peak spectrum can be described by a single power law with photon index V ~1.7 and has a flux of ~2 
photons cm -2 s _1 in the 20-200 keV energy band. The high sensitivity of IBIS has made it possible for the first 
time to perform detailed time-resolved spectroscopy of a GRB with a fluence of 7xl0 -6 erg cm -2 (20-200 keV). 
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1. Introduction 

Ever since their discovery ({Klcbesade l et al. 1973|l . 
Gamma-Ray Bursts (GRBs) have been a puzzling mys- 
tery, mostly because of their short durations and the ap- 
parent lack of counterparts at other wavelengths. A break- 
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through in this field came thanks to the Italian-Dutch 
satellite BeppoSAX, which had the capability to local- 
ize the bursts' prompt emission with a precision of a few 
arcminutes within a few hours. This led to the discov- 
ery of the afterglow emission at lower energies, initially 
in X-rays IjUosta et al. 1997|) and subsequently at opti- 
cal ( |van Paradijs et al. 1997| ) and radio IjFrail et al. 1997|) 
wavelengths, which allowed the redshift of these objects 
to be measured, and firmly established the cosmological 
nature of GRBs. 

Due to the limited duration and the fading char- 
acter of the afterglow emission, the prompt distribu- 
tion of GRB coordinates to the scientific community is 
a high priority. After the end of the BeppoSAK mis- 
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sion this task has been accomplished mainly by HETE- 
2 IjRicker et al. 2002jl . INTEGRAL IjWinkler et al. 1999|) . 
although not specifically designed as a GRB-oriented 
mission, can contribute to the rapid localization of the 
prompt emission of GRBs thanks to the INTEGRAL 
Burst Alert System (IBAS; |Mereghetti et al. 200ip . This 
software, running at the INTEGRAL Science Data Centre 
fISDC: ICourvoisier et al. 1999(1 . is able to detect and lo- 
calize GRBs with a precision of a few arcminutes in a few 
seconds, and to distribute their coordinates in near real 
time over the Internet. 

The high sensitivity of the INTEGRAL instruments 
also allows us to study in detail the prompt 7-ray emission 
of GRBs. This is particularly interesting for the faintest 
bursts, for which deep spectral studies were not possible 
up to now. For example, with the CGRO/BATSE instru- 
ment, time resolved spectroscopy was possible only for 
bursts with a fLuence larger than ^4xl0 -5 ergs cm~ 2 
IjPreece et al. 1998|) . 

On January 31 2003 at 07:38:49 UTC a GRB was de- 
tected in the field of view of the main instruments on 
board INTEGRAL: I BIS l|Ubertini et al. 1999jl and SPI 
IjVedrenne et al. 1999|) . Here we concentrate on the results 
obtained with IBIS, a coded mask imaging telescope based 
on two detectors, ISGRI and PICsIT, operating in the 15 
keV - 1 MeV and 170 keV - 10 MeV energy ranges, re- 
spectively. 

2. Detection and Localization 

GRB 030131 was discovered by IBAS (using IBIS/ISGRI 
data) on January 31 2003 at 07:39:10 UTC (~ 21 sec- 
onds after the beginning of the GRB, see below). The 
on-line automatic imaging analysis localized it to off-axis 
angles Z=10.1°, Y=-3.6°, in the partially coded field of 
view (only ~23% of the detector was illuminated by the 
GRB). The GRB coordinates were not distributed auto- 
matically by IBAS because most of the burst occurred 
during a satellite slew (IBAS is disabled during satel- 
lite slews). In fact only the first ^20 seconds of this 
~150 s long burst, during which the satellite attitude 
was stable and well known, were analyzed by IBAS, re- 
sulting in a low significance of the trigger. An off-line 
interactive analysis confirmed the reality of the event 
IjBorkowski et al. 2003|l . but the reported error region ra- 
dius was underestimated. A correct localization with an 
error radius of 5' (Mcrcghc tti et al. 2003a| ) was distributed 
only three days later. 

By accumulating data over short time intervals and an- 
alyzing the corresponding images, we confirmed that the 
satellite slew started at 07:39:09 UTC, as indicated by the 
attitude data. We therefore used the first 20 s of the event, 
corresponding to the stable pointing period, to derive the 
GRB position a J20 oo = I3 h 28 m 21 s , <5 J200 o = +30° 40' 
33", with an error radius of 2.5'. Although the statistical 
error in these coordinates is only 1.6', we conservatively 
added a systematic uncertainty of 2', based on the results 
obtained in IBIS observations of sources with known po- 
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Fig. 1. Localizations of GRB 030131: the annulus ob- 
tained by the IPN is consistent with the localization given 
in |Mereghetti et al. 2003a] (GCN 1847) and the one de- 
rived in this paper (the annulus obtained using SPI ACS 
and Ulysses is consistent with the one plotted but has a 
larger width). The cross indicates the position of the OT. 



sitions. Our final position for GRB 030131 is consistent 
with the one reported earlier fMcrcgh etti et al. 2003a) 
and with the annulus derived with the IPN using Ulysses 
and IBIS/ISGRI data (see Fig.[TJ. 

A provisional identification of an optical transient 
(OT) for GRB 030131 was reported by Fox et al. 
(2003a). The candidate OT, detected with the Palomar 
48-inch Oschin telescope + NEAT Camera, had magni- 
tude R^21.2 at 3.62 hours after the burst, but it was 
much fainter and barely detectable (R>23.5) ^26.8 hours 
after the burst at the 200-inch Hale telescope. Its coordi- 
nates, a J2 ooo = 13 h 28 m 22.29 s , <5 J20 oo = +30° 40' 23.7", 
are 20" from the center of the error circle derived here. 

As a follow-up, we obtained a 3x300 s exposure in the 
V band using the European Southern Observatory Very 
Large Telescope (VLT) with the FORS1 instrument at a 
mean date of 13 February 09:11:54 UTC. The seeing was 
about 1". There was no detectable object at the location 
of the candidate optical counterpart, with a 5 a upper 
limit of V>26.4. 

The marginal detection ~29 hours after the burst, 
with B = 25 .4±0.3 HGorosabel et al. 20031 but see also 
IHenden 20 03). and our VLT upper limit, confirm that this 
object is the OT of GRB 030131. Thus GRB 030131 is the 
first GRB detected with INTEGRAL with an associated 
optical counterpart. 
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GRB 030131 




Time after 07:38:49 UTC [seconds] 

Fig. 2. IBIS/ISGRI light curve of GRB 030131 in various 
energy bands ((a): 15-50 keV, (b): 50-300 keV, (c): 15-500 
keV). The six data gaps are artifacts caused by satellite 
telemetry saturation. Four peaks can be identified, after 
~15, ^40, ~55, ~85 seconds. The light curve has been 
corrected to take into account the varying fraction of the 
exposed instrument area during the satellite slew. Panel 
(d) shows the spectral variation of the GRB with time. 



3. Temporal and Spectral Analysis 

We have analyzed IBIS/ISGRI single events, for which 
arrival time, energy deposit and interaction pixel of the 
detector are known for each event. 

Fig.EJshows the light curves of GRB 030131 binned at 
1 s resolution in different energy bands. The burst started 
at 07:38:49 UTC and lasted for about 150 s. The time pro- 
file shows several peaks (note that the small gaps are ar- 
tifacts caused by satellite telemetry saturation). The Tg 
duration of the GRB in the 15-500 keV band is 124 s. 
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Fig. 3. Positions of GRB 030131 in detector coordinates 
as a function of time. Note the different scales of the two 



axes. 



Since the GRB peaks during the satellite slew, we 
could not use the instrumental coordinates to extract the 
peak spectrum. Therefore we made an image selecting a 
time interval of 1 s around the GRB main peak (t = 54 s 
in Fig- O - The high count rate at the peak allowed us to 
firmly establish the detector coordinates of the GRB even 
with this short integration time and thus to extract its 
peak spectrum. Since IBIS/ISGRI is a coded mask imag- 
ing instrument, the background can be estimated simulta- 
neously with the source flux, using the Pixel Illumination 
Function (PIF; Skinn er" 1995(1 . The spectra have been ex- 
tracted computing one PIF for each energy bin (128 lin- 
early spaced bins have been used between 19 keV and 1 
MeV) . Since a fully calibrated spectral response matrix for 
sources at large off-axis angles is not yet available, we di- 
vided the count spectrum by the closest (in detector coor- 
dinates) count spectrum of the Crab Nebula. The resulting 
photon spectrum can be well fitted by a power law model 
(x7 d °f = 8.66/8) with photon index T = l-73±g;if (90% 
confidence level). The flux is ~1.9 photons (~ 1.7xl0 -7 
erg) cm~ 2 s _1 (~6.5 crab) in the 20-200 keV energy range. 

To properly extract the total GRB flux and spectrum, 
we derived the GRB detector coordinates at various time 
intervals for the entire duration of the event. The first 
interval, corresponding to the stable pointing, lasts 20 s. 
The following 30 intervals last 3 s each and, finally, for 
the faint tail of the burst, four intervals with durations of 
5, 5, 10 and 20 s were used. The coordinates as a func- 
tion of time are shown in Fig. where one can see the 
apparent drift of -2° of the GRB in the IBIS field of 
view during the slew. To obtain the total (time averaged) 
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GRB 030131: Average Spectrum 



4. Discussion 
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Fig. 4. Time averaged IBIS/ISGRI spectrum of GRB 
030131. Data and best fit model are shown. 



spectrum, the 35 spectra extracted for the time intervals, 
corresponding to a net integration time of 147.2 s, were 
co-added. We derived the photon spectrum (shown in Fig. 
|IJ| in the same way that we did for the peak spectrum. The 
GRB is clearly detected up to 200 keV. In this case a sin- 
gle power law model does not provide a satisfactory fit 
(x 2 /dof = 39.84/18), indicating that a model which in- 
cludes a spectral break would give a better representation 
of the data. A fit using the Band model IjBand et al. 1993"|) 
yields a break energy .Erj=70±20 keV, a low-energy power 
law index a = 1.4±0.2 and a high-energy photon index 
/3=3.0±1.0 (x 2 /dof = 22.54/16). The fluence in the 20-200 
keV band is ~73.8 photons cm" 2 (7 x 10" 6 erg cm -2 ). 
To derive this value we have assumed that during the 4 
central telemetry gaps the GRB had its average spectrum 
and intensity. Extrapolating the spectrum to the BATSE 
energy range (50-300 keV) we obtain a fluence of ~23 pho- 
tons cm" 2 (3.2 x 10 -6 erg cm" 2 ). The values of peak flux 
and fluence derived above are consistent to within a factor 
of 2 with the ones measured with Ulysses. This indicates 
that our method does not suffer from large systematic er- 
rors and that the short telemetry gaps do not influence 
these values significantly. 

We have also investigated the spectral evolution. As 
in the case of the peak spectrum, the data for the indi- 
vidual time intervals can be fitted with a single power 
law, without evidence for a spectral break. The photon 
index as a function of time is plotted in Fig. A clear 
hard-to-soft evolution is seen. A hardening trend can also 
be seen corresponding to the rise of the second and third 
(main) peak with a softer spectrum for the latter peak. 
This correlation between light curve peaks and spectral 
hardening has already been reported in other bursts (e.g. 
IForri et al. 1993. 



The IBIS/ISGRI time-resolved spectroscopy of GRB 
030131 is consistent with the overall hard-to-soft evolution 
observed with BATSE in many brighter GRBs, for which 
this kind of analysis was possible (e.g. lPreece et al. "1 998 
IFord et al. 1995)l . The fluence of GRB 030131 is an or- 
der of magnitude smaller than those of the bursts stud- 
ied by those authors, indicating that such spectral be- 
haviour applies also to fainter GRBs. Clear evidence of this 
was also reported in the GRBs studied with BeppoSAX 
in the 2-700 keV energy range (e.s. IFrontera et aT~2 000 
IFrontera et al. 2003|) . While BATSE could better con- 
strain the break energy and the high-energy slope of the 
Band function, thanks to its higher relative sensitivity 
above 200 keV, IBIS/ISGRI allows more detailed stud- 
ies of the low-energy part of the spectrum for relatively 
low fluence GRBs. 

In the framework of the internal fireball shock 
model (Rces & Mcszaros 1994), and in particular of the 
Synchrotron Shock Model (Tavani 199(jJ), the hard-to-soft 
evolution can be interpreted in two ways. The first possi- 
bility is a decrease of the magnetic field in the postshock 
region as a consequence of the postshock flow expansion; 
the second is a postshock decrease of the index of the par- 
ticle distribution function as a consequence of strong cool- 
ing processes affecting the particle energy distribution for 
dynamical flow times larger than the radiating timescale. 
The two effects are not distinguishable in our case, since 
the time-resolved spectra do not have enough statistics to 
constrain a spectral break and hence a low-energy and a 
high-energy spectral index. The soft-to-hard evolution ob- 
served during the rise of the individual peaks, on the other 
hand, can be caused by an increase of the local magnetic 
field at the shock. Several authors reported that the dura- 
tion of single pulses in GRB time histories is energy depen- 
dent (e.g. |Link, Epstein fc Priedhorski 1993| ), with longer 
durations at lower energies, resulting in a hardening of the 
spectra before the peaks and a softening afterwards. 

The optical transient associated with GRB 030131 in- 
dicates that we can classify it as an "optically dim" GRB. 
In fact it is as faint as (or even fainter than) the tran- 
sient associated with GRB 030227, which is the only other 
INTEGRAL GRB with a firmly established optical con- 
terpart, and was detected at R~23.3 12 hours after the 
burst | |Mereghetti et al. 20 03b). It is also comparable to 
GRB 021211. This event is also considered an optically 
dim burst since it was detected at R~18.2 1.3 hours after 
the prompt emission IjFox et al. 2003 b|l and was fainter 
than R-22.5 after 12 hours <|Klose et al. 2002jl . In ad- 
dition GRB 030131 is located at much higher Galactic 
latitude (b ~ 81°) which implies smaller foreground op- 
tical extinction. This indicates that, despite the efforts 
of observers, in some cases, optical follow-up with small 
telescopes is not an easy task, even less than 1 day af- 
ter the burst (e.g |Fynbo et al. 2001] ). The prompt lo- 
calization of GRBs is hence a high priority in order to 
achieve a successful follow-up. The results on GRB 030501 
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ijEeck" mann et al. 2003|) (an alert 30 s after the start of the 
GRB with an uncertainty of 4.4'; IBAS Alert 596) show 
that IBAS is now able to provide this service. 
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